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Abstract 

Background: Assessment of dual time point (DTP) positron emission tomography was carried out with the aim of a 
quantitative determination of K m , the metabolic uptake rate of [ 18 F]fluorodeoxyglucose as a measure of glucose 
consumption. 

Methods: Starting from the Patlak equation, it is shown that K m & m t /c^ + V r /r ai where m t is the secant slope of the 
tissue response function between the dual time point measurements centered at t = to. c° Q = c a (to) denotes arterial 
tracer concentration, V r is an estimate of the Patlak intercept, and z a is the time constant of the c a (t) decrease. We 
compared the theoretical predictions with the observed relation between K s = rnt/c® and K m in a group of nine 
patients with liver metastases of colorectal cancer for which dynamic scans were available, and K m was derived from 
conventional Patlak analysis. Twenty-two lesion regions of interest (ROIs) were evaluated. c a (t) was determined from a 
three-dimensional ROI in the aorta. Furthermore, the correlation between K m and late standard uptake value (SUV) as 
well as retention index was investigated. Additionally, feasibility of the approach was demonstrated in a whole-body 
investigation. 

Results: Patlak analysis yielded a mean V r of V r = 0.53 ± 0.08 ml/ml. The patient averaged x Q was 99 ± 23 min. Linear 
regression between Patlak-derived K m and DTP-derived K s according to K s = b • K m + a yielded b = 0.98 ± 0.05 and 
a = —0.0054 ± 0.001 3 ml/min/ml (r = 0.98) in full accordance with the theoretical predictions b = 1 and a & —V r /r a . 
K s exhibits better correlation with K m than late SUV and retention index, respectively. /C s (c) = K s + V r /r a is proposed as 
a quantitative estimator of K m which is independent of patient weight, scan time, and scanner calibration. 

Conclusion: Quantification of K m from dual time point measurements compatible with clinical routine is feasible. 
The proposed approach eliminates the issues of static SUV and conventional DTP imaging regarding influence of 
chosen scanning times and inter-study variability of the input function. K s and exhibit improved stability and 
better correlation with the true K m . These properties might prove especially relevant in the context of radiation 
treatment planning and therapy response control. 
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Background 

For many years, quantification of the metabolic rate of glu- 
cose consumption with dynamic [ 18 F]fluorodeoxyglucose 
(FDG) positron emission tomography (PET) using the so- 
called Patlak plot, a procedure most clearly described by 
Patlak in his seminal papers [1,2], has proven valuable in 
PET research and clinical routine. 

However, in the clinical oncological setting, quantifica- 
tion is mostly restricted to the ubiquitously used standard 
uptake value (SUV). The reason is twofold: (1) no need 
(or even inability) to determine the arterial input function 
(AIF) and (2) inability to perform dynamic whole-body 
investigations. 

Without question, the SUV (defined as the tracer uptake 
at a certain time point normalized to injected dose per 
unit body weight) has proven a valuable means of achiev- 
ing a certain level of quantitative description, thus allow- 
ing, e.g., definition of standardized evaluation schemes 
(see [3] for an overview). 

The approach, however, has known shortcomings [4-6] . 
SUVs do not directly provide information about the tracer 
kinetics but, by their very nature, only a static snapshot 
somewhere on the tissue response function (TRF). Natu- 
rally, SUVs are varying along the given TRF and are thus 
prone to variability when not determined at a strictly stan- 
dardized time. Since SUVs do not contain any information 
of the actual rate of tracer accumulation (related to the 
slope of the tissue response function), TRFs from different 
tissues might in extreme cases even intersect at a certain 
time (thus exhibiting identical SUVs and zero image con- 
trast at this moment) while having completely different 
kinetic properties. Tissue SUV stability is further com- 
promised by not accounting for the sizable inter-study 
variability of arterial blood SUV which directly influences 
the actually obtained tissue uptake. 

One quite extensively investigated way around the snap- 
shot problem' is dual time point (DTP) investigations [7,8] 
in which two successive whole-body scans are performed 
to obtain information regarding the rate of tracer accumu- 
lation. While being undoubtly valuable in discriminating 
between tumor and inflammation, quantitative evaluation 
of DTP measurements is usually restricted to computa- 
tion of a so-called retention index, RI, representing the 
percentage change of SUV max or SUV mea n between early 
and late images (see, e.g., [9,10]). However, the retention 
index, too, depends on the acquisition time of (and time 
difference between) early and late PET scan and, there- 
fore, requires the same strict standardization as the SUV 
approach to provide useful quantitative measures. The 
retention index, too, is affected by the mentioned AIF 
variability at late times due to the evoked changes of the 
TRF slope. 

There also have been attempts to directly use the TRF 
slope obtained in dynamic scans as a substitute for actual 



kinetic modeling [11,12] while avoiding measurements 
of tracer concentration in blood. However, a convinc- 
ing physiological interpretation of the slope parameter 
is missing. Furthermore, the approach suffers from the 
same problems as SUVs and retention index regarding 
the uncontrolled influence of the inter- and intra-subject 
variability of the AIF. 

In this study, we propose a new assessment of DTP (and, 
more generally, TRF slope) -based methods with the aim of 
a quantitative determination of K m , the metabolic uptake 
rate of FDG. We demonstrate that starting from the Pat- 
lak model, one can derive an analytical relation between 
K m and the TRF slope m t , which only requires the image- 
based determination of the AIF during the respective late 
PET scans. The derived relation is especially compatible 
with dual time point whole-body investigations. 

In this retrospective investigation, we evaluate the new 
approach in a group of patients with liver metastases of 
colorectal cancer for which K m was determined, both, 
by conventional Patlak analysis of the fully dynamic PET 
scans as well as by the newly developed approach. 

Methods 

Theory 

It is well known that the TRF after a bolus injec- 
tion of FDG appears to be approximately linear at later 
times. Closer inspection, however, reveals, that the curve 
exhibits a finite curvature: the slope decreases with time 
due to the continuously decreasing AIF (see Figure 1). 
In the Appendix, we demonstrate that for times t when 
the Patlak equation is valid (usually for t > 20 — 30 
min), the ratio between the instantaneous values of TRF 
slope and AIF level can be expressed in terms of the 
parameters K m and V r of the Patlak model and the time 
constant x a describing the essentially mono-exponential 
decrease of the AIF in the considered time window. It 
is shown in the Appendix that the TRF slope at t = to 
is very nearly identical to the slope of the secant con- 
necting the boundary points of a finite symmetric time 
interval around to (and also to the average slope in this 
interval). 
One finally arrives at the relation 

K m =K s + ^ = ^ + ^ (1) 

ta c a t& 

with 



where m t is the secant (or average) TRF slope in the cho- 
sen time interval centered at to and <P a = c a (fo) (see 
Figure 1). 

The rate K s defined by Equation 2 (i.e., the ratio between 
the TRF slope and AIF level at time to) can be deter- 
mined from measurements during the late phase alone. 
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Figure 1 AIF plus TRF calculated for K\ = 0.3 ml/m in/ml, ki = 0.5 /min, A3 = 0.08 /min. The TRF does not become linear at later times but 
exhibits a visible curvature. However, the slope at some time point to (to = 60 min in this example) is nearly identical to the slope of the secant 
connecting the boundary points of a finite time interval centered at t 0 .The AIF is scaled such that c a (to) = I.The data accessible in a DTP 
measurement are indicated by the square plotting symbols. For further details, see the main text. 



Contrary to the Patlak method, knowledge of the full AIF 
is not required. To the extent that K s ^> V r /x a) K s might 
directly serve as an (negatively biased) approximation of 
K m . Moreover, to the extent that V r can be replaced by a 
suitable constant value V r , K s differs from K m only by a 
x a - dependent offset that can be added to K s to obtain a 
corrected value 

=K S +K 0 =K S + — (3) 

la 

that approximates K m quite accurately (see Appendix and 
Figure 1). 

We have compared these theoretical predictions with 
the actually observed relation between K s and K m in a 
group of patients with liver metastases for which fully 
dynamic scans were performed. 

Study sample 

The investigated patient group included nine male sub- 
jects with liver metastases of colorectal cancer (mean 
age 62.8 years, range 48 to 76). For each patient, one to 
three dynamic PET scans of 60 min duration were per- 
formed (altogether 15 scans). Scans started immediately 
after injection of 346 to 430 MBq FDG. The scans were 
performed with an ECAT EXACT HR+ (Siemens/CTI, 
Knoxville, TN, USA). The acquired data were sorted into 
23 to 31 frames with 10 to 20 s duration during bolus 
passage, 30 to 150 s duration until 10 min post-injection 
(p.L), and 300 s duration afterwards. Tomographic images 
were reconstructed using attenuation-weighted OSEM 
reconstruction (6 iterations, 16 subsets, 6 mm FWHM 
Gaussian filter). 



Additionally, feasibility of the generation of paramet- 

(c) 

ric K s maps was demonstrated in a whole-body FDG 
investigation of a 63-year-old woman with bronchial car- 
cinoma of the right lung and lymph node metastases of the 
right hilar region and the mediastinum (Philips Ingenuity 
TF PET/MR (Philips, Cleveland, OH, USA), injected dose 
273 MBq, first scan 67 min p.i. (2 min per bed position), 
second scan 117 min p.i. (1.5 min per bed position)). 

Data evaluation 

Region of interest (ROI) definition was performed using 
ROVER (ABX, Radeberg, Germany) [13,14]. The AIF was 
determined from a roughly cylindrical three-dimensional 
(3D) ROI centered in the aorta using a concentric safety 
margin of at least 1 cm to exclude partial volume effects. 
3D lesion ROIs were denned in 22 lesions, and the respec- 
tive TRFs were computed. Further data analysis was 
performed using the R software for statistical computa- 
tion [15]. 

For all 22 lesions, K m and V r were derived from the con- 
ventional Patlak analysis of the full dynamic data later than 
20 min p.i. (at which time all Patlak plots already were 
linear). For comparison with the corresponding result of 
the subsequent DTP evaluation, x a was determined from a 
mono-exponential fit to the complete AIF data in the time 
window used for the Patlak analysis. Variability of x a and 
V r was expressed as mean ± standard deviation (SD). 

Dual time point data were generated from the data 20 
to 30 min and 50 to 60 min p.i., yielding two pairs of c a 
and c t values which were assigned to the respective frame 
centers = 25/55 min (which corresponds to to = 

0.5 • (t~ + f+) = 40 min and At = t+ - t~ = 30 min). 
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Using the abbreviations Cy a = Ct/ai^) and Equation 2, K s 
is given by 



- - - 



Ac t 
~At 



(4) 



where c Q a was calculated from the exponential connect- 
ing the two points (t~, c~) and c+) which yields C® = 

For computation according to Equation 3, we fixed 
V r to the mean of the Patlak evaluation for all lesions 
(V r = 0.53 ml/ml), while x a was estimated individually for 
each study from the exponential connecting (t~,c~) and 
(£+c+)as 

At 
In (Ca/ct) 

Additionally, the retention index was computed as 
RI = Act/Cj". Linear regression analysis was performed 
between K m and K s , K^ c \ RI, and c+ (the SUV of the 
lesions in the late image), respectively. Parametric images 

(c) 

of K s ,I<r, and/C m were generated for visual comparison 
after filtering of the DTP image data with a bilateral fil- 
ter [16] (spatial filter width 9 mm, intensity filter width 2.5 
SUV). 

Influence of image noise 

Considering a single voxel and neglecting the (much 
smaller) statistical error of the ROI-based c° a value, it fol- 
lows from Equation 4 that the relative statistical error of 
K s is equal to that of Ac t and thus, by Gaussian error 
propagation, 



M + M 



which decreases with increasing concentration difference 
Ac t . Taking into account that measurement times of both 
dual time point measurements might be adjusted in such 
a way that a c - ~ o c + , one can get a rough estimate of the 
error according to 

1 G ct 



V2- 



1 



where the final ratio represents the relative SUV error of 
the second dual time point measurement. For At & 30 
min and typical tumor accumulation rates of & 2% to 4% 
per minute, one can thus estimate that the relative errors 
of K s are about 2.5 to 4 times higher than the correspond- 
ing SUV errors (the statistical error of is quite similar 
since the small correction term V r /x a cannot contribute 

(c) 

much to the total statistical uncertainty of IQ ). Although 
the noise in the parametric maps can thus be expected to 
be distinctly higher than that in the uptake images, the 



resulting visual quality is still quite satisfactory for reason- 
able choices of At (> 30 min) as will be demonstrated in 
the following. 

Results 

The obtained results are summarized in Tables 1 and 2. 
Figure 2A shows the correlation between Patlak-derived 
K m and K s . The solid line is the line of identity, and the 
dashed line is the linear regression result. The linear cor- 
relation is very good, and the fitted slope is identical to one 
within the given error limits of about 5% (Table 2). The 
fitted intercept of —0.54 ml/min/100 ml thus represents 
the experimentally observed average underestimate of the 
true K m by K s . 
Figure 2B presents the correlation between K m and 

(c) 

IQ = K s + Ko according to Equation 3. The correction 
term K$ = V r /r a was computed using the average V r 
derived from the Patlak analysis of all 22 lesions, V r = 0.53 
ml/ml, and individual (investigation-specific) time con- 
stants x a derived from the DTP data (the independent 
determination of x a from the full dynamic data in the 
Patlak time window - performed as a consistency check - 
yielded essentially the same result (104±20 min (dynamic) 
vs. 99 =b 23 min (DTP)) but was not used further). As 
can be seen, the degree of linear correlation is distinctly 
improved in comparison to Figure 2A. Furthermore, the 
fitted straight line now essentially coincides with the line 
of identity. Consequently, the average difference between 

(c) 

K m and IQ amounts to only 1.4 ± 4.1% and exceeds 10% 
only in a single lesion. 

For comparison, Figure 3A,B presents the correlations 
between K m and the late SUV uptake cf, and K m and the 
retention index RI, respectively. Obviously, the correla- 
tion between K m and c~l is rather poor. The correlation 
between K m and the RI is substantially higher but still 
clearly below the degree of correlation between K m and K s 
or Ks C \ The correlation in Figure 3A is clearly distorted by 
the group of the six highest observed SUV values which 
correspond to only moderately high K m values. This phe- 
nomenon might be explained by the exceptionally high 
c° a values observed in the respective patients (see inset 
graphic in Figure 3A). 

Table 1 Summary of parameters entering the K s and ff s (c) 
determination 



Mean ± SD 



Range 



V r (ml/ml) 


0.53 ± 0.08 


0.39- 


-0.68 


Dynamic r a (min) 


1 04 ± 20 


79- 


156 


DTP x a (min) 


99 ±23 


81 - 


172 


c° a (SUV) 


3.1 ±0.7 


2.4- 


■4.5 



'Dynamic r a ' denotes the result derived from the full dynamic AIF data with 

r > 20 min. 'DTP r a ' was computed solely from the used DTP data. Only the latter 

was used in further computations. 
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Table 2 Linear regression results: Pearson correlation 
coefficient r and the obtained regression parameters are 
shown 





r 


Slope (mean ± SD) 


Intercept (mean ± SD) 


Km vs. K s 


0.98 


0.98 ± 0.05 


-0.54 ±0.1 3 




0.99 


0.98 ± 0.04 


0.04 ±0.10 


Km VS. C+ 


0.60 


1 .62 ± 0.48 


2.81 ± 1.33 


K m vs. Rl 


0.88 


1 1 .38 db 1.34 


5.44 ±3.79 



All correlations were found to be statistically significant with p < 0.005. Chosen 
units are in milliliters per minute per 100 ml for K m ,K s , and K^; SUV for c+; and 
percentage for Rl (from which the respective slope and intercept units follow). 



Figure 4A provides one example of a lesion uptake 
image, and the corresponding parametric images of K s , 

(c) 

IQ , and K m are shown in Figure 4B,C,D, respectively, all 
of which are displayed in a common scale. The compara- 
ble, enhanced target-to-background contrast of the three 
parametric images relative to the uptake image is obvious. 
In agreement with the theoretical expectation and the ROI 
data in Figure 2, there is good quantitative concordance 
between K m and Ks°\ while K s exhibits a constant negative 
bias of about Kq ^ —0.5 ml/min/100 ml in comparison to 
K m . 

Finally, Figure 5 demonstrates the feasibility of generat- 

(c) 

ing parametric IQ maps of reasonable statistical quality 
for a typical DTP whole-body FDG study. 

Discussion 

Our main result is that in the investigated patient group, 
there is a very pronounced linear correlation K s = a + b • 
K mi where b is very nearly equal to one (see Figure 2A). 
This behavior is in complete agreement with the formal- 
ism presented in the Appendix, notably Equation 13: the 
variations of the (small) term V r / x a should be essentially 
uncorrected to K m so that a high (but slightly "noisy") lin- 
ear correlation between K m and K s with a slope near one 



is predicted. Furthermore, according to Equation 13, the 
modulus of the intercept, a = —0.0054 ml/min/ml, should 
be approximately equal to the average of V r /r a in the 
investigated patient group. This prediction, too, is in com- 
plete agreement with the actual values of V r (determined 
from Patlak analysis) and x ai namely V r = (0.53 ± 0.08) 
ml/ml and x a = (99 ± 23) min. 

The second important finding is the fact that the degree 
of correlation as well as quantitative agreement between 
K m and K s can be further improved by assuming a rea- 
sonable constant value for V r (since V r is inaccessible 
in DTP measurements) and determining individually the 
rate of decrease, x a , of the AIF (which can be esti- 
mated from the DTP measurement). The corrected K s , 

(c) ~ 

IQ = K s + V r /r a , exhibits an improved correlation 
to K m (due to compensation of the x a variability) and 
also improved quantitative concordance as long as V r is 
roughly in accord with the individual true V r . Both phe- 
nomena are illustrated clearly in Figure 2B. The distinctly 
improved correlation (compared to Figure 2A) is achieved 
by the individual correction of the x a influence. The resid- 
ual deviations from the perfect correlation in Figure 2B 
are mainly due to the variability of V r . A nearly perfect 
quantitative agreement with K m is observed since V r was 
set to the mean of the actual V r values derived from 
the Patlak analysis. This obviously would not be possible 
when considering realistic DTP measurements (without 
a preceding complete dynamic study), and a less-than- 
perfect quantitative agreement should be expected in this 
case. Nevertheless, as the comparison of Figure 2A,B 
suggests, performing the correction with some roughly 
correct value for V r will always decrease the bias between 
K s and K m . 

The rather small variability of tumor V r observed in the 
present investigation might seem surprising. However, the 
square of ^2/(^2 + £3) appearing in Equation 6 will never 
deviate very much from unity since for FDG, I<2 quite 




K m (ml/min/1 00ml) K m (ml/min/1 00ml) 



Figure 2 Correlation between the metabolic rateK m and the DTP-derived rate constants. (A) K m and K s and (B) K m and Kf } . Solid black lines 
represent the line of identity; dashed lines represent the least squares straight line fits to the data. 
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Figure 3 Correlation between the metabolic xa\eK m and Standard Uptake Value and retention index, respectively. (A) K m and cf (in SUV 

units) and (B) K m and retention index (in percent), cf is the tissue uptake at f = 55 min. Dashed lines represent the least squares straight line fits to 
the data. The inset graph in (A) additionally provides the correlation between c° a and cf , which helps to explain the group of data points above the 
regression line with SUV>8 (see the 'Discussion' section). 



generally is distinctly larger than £3. The variability of V r 
is thus mostly controlled by the first term, K\/k2. Since 
both K\ and I<2 are usually identified as being associated 
with the facilitated diffusion across the cell membrane, 
it might very well be expected that the ratio I<i/I<2 is 
essentially constant, independent of the actual K\. This 
might be the underlying reason for the low variability 
of V r observed in this study. Whether V r variability is 
higher in other tumors remains to be investigated, but 
we believe this to be unlikely. V r should never be much 
larger than about 0.6 to 0.7 ml/ml which appears to 
be a rough upper bound for the Ki/I<2 ratio. Accord- 
ing to our own data, this is true, e.g., in the human 
brain (I<i/k 2 « 0.1/0.15 = 0.67 ml/ml) as well as the 
myocardium (I<i/k 2 « 0.6/1.4 = 0.43 ml/ml). V r in these 
organs is rather low 0.3 ml/ml) due to the large I<3 in 
both tissues. 

We surmise, therefore, that V r in tumors (and healthy 
tissue) will never deviate too much from the value of 0.53, 
ml/ml used for K s correction in this study. The corrected 

(c) 

rate, IQ , can then be expected to be a less biased estimate 
of K m than K s over a substantial range of actually realized 
V r values between about 0.2 and 0.7 ml/ml (see Figures 6 



(c) 

and 7). Whether IQ does offer any advantages over K s 
in terms of clinical relevance remains to be seen, but the 
improved correlation with K m seems justification enough 
to perform the correction. 

(c) 

The very high correlation between Ks (or K s ) and K m 
is to be compared with the markedly inferior correla- 
tion between K m and late SUV (cf) and retention index 
RI, respectively (Figure 3). Since all these parameters are 
ultimately intended as surrogate parameters of K m , the 
superiority of K s seems obvious. Since in the present study 
the retention index is computed from exactly the same 
DTP tissue data as K s , it is worth to point out that the sole 
factor responsible for the much better K s (K m ) correlation 
is adequate consideration of the substantial inter-subject 
c® a variability (see Table 1). Indeed, one could write K s = 
Ac t /At/c° a = Ac t /Ct/At • c-/c° a = RI/At • c~ /c% where 
At is just a constant in the present context. To some 
extent, K s might thus be considered just a more sensible 
definition of a retention index where the uptake difference 
Act is normalized to c° a (as well as A£) instead of cj . 

The observed very low correlation between late SUV 
and K m is caused by the six data points with SUV>8 
in Figure 3A. Leaving these six points out increases 




D 



Figure 4 Liver metastasis of a colorectal carcinoma exhibiting a central necrosis. (A) A representative sagittal slice of uptake and the 
corresponding parametric images of (B) K Sl (C) kP , and (D) K m , respectively, are shown. The means of the rate constants over the lesion are as 
follows: K s = 2.23 ml/min/1 00 ml, K$ = 2.80 ml/min/1 00 ml, and K m = 2.90 ml/min/1 00 ml, demonstrating a quite satisfactory quantitative 



agreement between and K, 
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Figure 5 Comparison of parametric images of SUV and 7<^ c) . (A) FDG SUV 67 min. p.i. and (B) parametric A* map in a 63-year-old woman. The 
patient has a bronchial carcinoma of the right lung and lymph node metastases of the right hilar region and the mediastinum (coronal maximum 
projection). Note the different regional contrast and the used units in both images. 



the correlation coefficient to 0.94 which is in good 
agreement with published data [17]. Closer inspection 
revealed exceptionally high c Q a values (see inset graphic 
in Figure 3A) for the affected data points which might 
have physiological reasons but could also hint at erro- 
neous SUV calibration (for which, however, a retrospec- 
tive inspection did not find any evidence). In any case, 
the data demonstrate the high sensitivity of SUV evalu- 
ations to variations of the AIF level and incorrect SUV 
calibration. 

The comparison of uptake and parametric images in 

(c\ 

Figure 4 demonstrates that K s as well as IQ reproduces 
the essential features of the Patlak K m image, notably 



the increased contrast between metastasis and liver back- 

(c) 

ground. Regarding the targeted lesions, the IQ ; image is, 
moreover, in good quantitative agreement with the K m 
image and could thus serve as a basis for regional quan- 
titative evaluation. We, therefore, believe it is worthwhile 

(c\ 

to investigate the potential suitability of K s as a quan- 
titative estimator (and not just a surrogate) of K m more 
thoroughly in future studies. Figure 5 demonstrates that 

(c) 

reasonable statistical quality of the IQ map can in fact be 
achieved in whole-body DTP investigations as well. 

Compared to more conventional approaches, our 
approach has several relevant benefits. The most impor- 
tant one in our view is the potential to perform fully 




Figure 6 Visualization of the difference between K m and K5. The plots cover a substantial range of the parameters kj and k^, assuming a fixed 
ratio K] /k 2 = 0.7 ml/ml (a rationale for fixing this ratio is given in the discussion above). Top left, K m ; top right, V r ; bottom left, absolute difference 
(Km - K s ); and bottom right, fractional difference {(K m - K s )/K m ). Parameters and their respective units: k 2l k 3 (1/min); V r (ml/ml); and K m , K s 
(ml/min/1 00 ml). Moving along the line, V r ~ 0.55 ml/ml between K m = 1 and 4 ml/min/1 00 ml corresponds approximately to the experimental 
data of Figure 2A. 
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Figure 7 Visualization of the difference between K m and Kg for ~ta = 99 min and an assumed distribution volume V r = 0.53 ml/ml. In 

comparison with Figure 6, the systematic difference relative to the true K m is mostly removed even when V r deviates distinctly from the assumed 
value. Moving along the line, V r ~ 0.55 ml/ml between K m = 1 and 4 ml/min/1 00ml corresponds approximately to the experimental data of 
Figure 2B. 



quantitative whole-body investigations based on a DTP 
acquisition. The only additional prerequisite is identifica- 
tion of the aorta or left ventricle in the DTP data. One 
gains the ability to directly identify regions of elevated 
irreversible FDG metabolism and to put the established 
DTP approach on a quantitative basis. A further advan- 
tage is the implied correction for the sizable inter-subject 
variation of the blood tracer concentration (SUV range, 
2.4 — 4.5 in this study). The latter correction alone clearly 
improves the correlation between the derived parame- 
ter (K s ) and the targeted one (K m ). Another important 
aspect is elimination of the dependence of SUV uptake 
and retention index on the time of measurement (s). To 
the extent that the Patlak model can be considered valid 
(negligible fcj), the proposed procedure yields a time- 
independent result, namely a direct estimate of the invari- 
ant rate K m which prospectively should allow definition 
of improved, objective reference values. A further impli- 
cation is elimination of any intra-scan time dependence 
in whole-body/multi-bed studies. Last but not least, the 
issue of ensuring correct SUV calibration is eliminated 
since all calibration factors cancel out when perform- 
ing an image-based determination of both TRF slope 
and c Q a . This observation seems especially relevant for 
multi-center studies. 



Conclusion 

We have demonstrated that it is possible to derive a 
quantitative estimate of K mt the metabolic trapping rate 
of FDG, solely from a dual time point measurement. 
We believe this approach to be of potential relevance 
especially in the context of oncological whole-body inves- 
tigations where the required AIF information is available 
in the field of view (aorta or left ventricle). In this case, the 
approach eliminates most if not all issues of static SUV 
and conventional dual time point imaging regarding the 
influence of the chosen scan times relative to the time of 
injection and the substantial influence of inter-study vari- 
ability of the AIF. Consequently, the derived parameters 

(c) 

K s and IQ exhibit a much improved stability and much 
better correlation with the true K m . These properties 
might prove especially relevant in the context of radia- 
tion treatment planning and therapy response control. 
Whether this is indeed the case has to be investigated in 
appropriate future studies. 

Appendix 

We start with the standard Patlak formula but avoid divi- 
sion by c a (t):* 

C t (t)=K m - f C a (s)ds+Vr-C a (t), (5) 

Jo 
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where K m is the metabolic trapping rate, defined by 
K - Klh 

and V r is the apparent volume of distribution defined by 



V r = 



K x k 2 Ki 



(I<2+I<3) 2 l<2 \k 2 +k3 



k 2 



(6) 



Equation 5 is valid for times t > T* where T* « 20 to 30 
min pi.. Utilization of this equation for K m determination 
requires measurements of the TRF only for t > T* but 
measurement of the complete AIF starting at time zero. 
We now want to eliminate the dependency on measure- 
ments prior to T*. By taking the time derivative at some 
time point t > T*, it follows directly from Equation 5 that 

c t (t) = K m • c a (t) + Vr-Ca(t) 

or after division by c a (t) (suppressing the t argument) 



= K m + V r 



Ca 



(7) 



Focusing on some specific time point t = to, we use the 
Taylor expansion of c a {t) around to {c^ (to): nth derivative 
at t = t 0 ): 



.w = E 



«=0 



(t - t 0 ) n . 



(8) 



Introducing the parameters x n denned by (— x n ) n 
c a/ c a^(to)> Equation 8 can be rewritten as 



(9) 



where ro is always equal to one. The parameters x n> o are 
constructed in such a way that for a mono-exponential 
decrease of c a (t) near to, we obtain x n> o = x a , where x a is 
the time constant of the exponential. Actually, it is known 
that starting rather early after bolus injection (t > 20 
min), c a (t) can be reasonably well described by a slow 
mono-exponential decrease with a time constant x a ~ 100 
min (in the present study, we found an average value of 
x a = 99 min, while a value of x a — 80 min was reported 
in [18]). 

Inserting the Taylor expansion from Equation 9 into 
Equation 7, we get (c® = c t (to)) 



t 



(10) 



In order to derive K m from this equation, we need to 
reliably estimate c®/c® as well as to have knowledge of 



1/ti (the fractional rate of decrease of the AIF at t = 
to). Obviously, direct determination of the time deriva- 
tive c® at t = to is not feasible in real (noisy) data. On 
the other hand, it is not clear whether the average slope 
over a necessarily rather large neighborhood (required for 
reasons of limited time resolution and count rate statis- 
tics) is an acceptable approximation of c® (since the slope 
changes over time). For investigation of this question, 
we compute from Equation 5 the difference c\ — = 
ct (to + ^) — ct (to — ^p) for two time points lying sym- 
metrically around to at a finite (possibly large) distance At 



Act(At) =cf-c t =K m - c a (s)ds+ V r • (c+ - c a ) 

(11) 

with c± = c a (t 0 ±f). 

Replacing all occurrences of c a (t) in Equation 11 
by the Taylor series in Equation 9 (neglecting fourth 
and higher order terms) and executing the integra- 
tion separately for each term of the series yield 
after some straightforward but lengthy calculations the 
following equation: 



Act = K, 



1+ 24 \x^) 



Vr 



1 + 



111 

24 r 3 



(12) 



At 



The detailed derivation of Equation 12 is presented in 
an additional file (see Additional file 1). The factors in 
square brackets deviate only minimally from one up to 
even quite large values of At. For the sake of simplicity, we 
will demonstrate this only for the well-established approx- 
imately mono-exponential behavior of c a (t) at later times 
but emphasize that the conclusions remain the same when 
using other reasonable parametrizations of the observed 
shape of the AIF at later times (e.g., by an inverse power 
law). 

As already pointed out, for a mono-exponential 
decrease of c a (t), all x n> o coincide with the time constant 
x a of the exponential. Consider, then, choosing At — 60 
min in Equation 12 . Since x a ^ 100 min, we have for both 
square brackets 1 + 1/24 • 0.6 2 = 1.015. It is, therefore, 
permissible to replace both square brackets by one. This 
yields 



Ac t 



Vr 



c" • At. 



Thus, Ac t is to a very good approximation proportional 
to At. At can become quite large, e.g., At = 1 h, as long as 
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the lower bound to — 4^ remains larger than T*. Defining 
the secant slope m t between and 



m t = 



Ac t 
At 



and introducing the rate constant K s 

for the ratio of the secant slope and the blood concentra- 
tion at to, we get 



V, 



mt 

K s = -^= K, 



K m =K S + 



Vr 



(13) 



(14) 



Comparison of Equation 13 with Equation 10 yields the 
important result 

m t = c° t . 

In other words, the secant slope is to a very good approx- 
imation equal to the instantaneous slope at to and thus 
can be used instead. This in turn implies that the aver- 
age slope of the TRF (derivable, e.g., by a least squares fit 
of a straight line in the considered time window), too, is 
very nearly identical to m t . Note that these conclusions are 
valid even if c t (t) varies considerably over the considered 
time interval (see Figure 1). Formally, this result is iden- 
tical to stating that a second-order Taylor expansion of 
c t (t) around to turns out to be sufficiently accurate within 
to ± At/2. 

The quantitative relation between K s and K m is inves- 
tigated in Figure 6. For this figure, we computed K m and 
V r over a range of sensible choices for the transport con- 
stants K\, I<2, and £3. The resulting K m and V r (top row of 
Figure 6) are used to compute K s from Equation 13 for a 
realistic value of x a (we chose x a = 99 min). The bottom 
row in Figure 6 compares the true K m to K s . 

As can be seen (bottom right), the fractional devia- 
tion of K s from K m becomes large only when I<3 is very 
small (i.e., when there is virtually no trapping). Overall 
K s is a negatively biased estimator of K mt but an approx- 
imate correction of the bias is possible considering the 
following. 

According to Equation 14, theVr and K m — K s maps 
in Figure 6 differ only by a constant factor x a (and a 
conversion factor of 100 due to the chosen units of 
ml/min/lOOml for K m and K s ). Moreover, V r does vary 
only modestly in comparison to the individual rate con- 
stants and to K m (except when £3 becomes distinctly larger 
than I<2, but this is not observed in real data). Therefore, 
K m — K s does not vary much across the relevant part of 
the /<2//<3 plane. We, therefore, hypothesize that the differ- 
ence K m — K s can be actually treated to be approximately 



constant. Consequently, we propose to estimate K m using 
only late time measurements of c a (t) and c t (t) as follows: 

1. Determine the secant TRF slope m t in the time 
interval t^ = to =b ^ from a dual time point 
measurement of c t (t) starting at sufficiently late 
times after injection, typically t > (20 — 30) min. 

2. Estimate c Q a = c a (to) and r a from the exponential 
connecting the two time points t~ ,t + . 

3. Compute K s = m t /c^. 

4. Compute a correction term Kq = V r /r a using the 
individually determined x a and a fixed value V r for 
the distribution volume. In the absence of any 
specific information regarding V r in the investigated 
tumor entity, we propose to use the average V r 
determined in this study, i.e., V r = 0.53 ml/ml. 

5. Finally, compute the corrected K s , i.e., 



/C 5 (C) =K s +Ko=Ks + V r lx a 

as a quantitative estimate of the true K n 



(15) 



(c) 

According to Equations 14 and 15, Ks is equal to K m if 
V r = V r (irrespective of the values of K±- !<% yielding this 

(c) 

V r value). Therefore, K s remains a very good approxima- 
tion of K m as long as V r does not deviate too much from 
the assumed value. This behavior is illustrated in Figure 7. 

Endnote 

a For completeness, we mention that in the presence of 
sizable fractional blood volume (fbv), the substitutions 
K m -> K m = (1 -fbv)K m and V r -> V? = (1 -fbv)V r + 
fbv = V r + fbv(l — V r ) would have to be performed in 
Equation 5 where the asterisked' quantities would be the 
experimentally accessible ones. 

Additional file 



Additional file 1: Derivation of Equation 12. A pdf file showing the 
complete derivation of Equation 1 2 using Taylor expansion. 
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